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Abstract 

The BL Lac Object Mkn 501 was in a state of high activity in the TeV range in 1997. During this time Mkn 
501 was observed by all Cherenkov-Telescopes of the HEGRA-Collaboration. Part of the data were also taken 
during moonshine thus providing a nearly continuous coverage for this object in the TeV-range. We have 
", carried out a QPO analysis and found evidence for a 23 day periodicity. We applied the same analysis on 

Q\ ' the 'data by dwell' x-ray lightcurve from the RXTE/ASM database and found also evidence for the 23 day 

^ ■ periodicity. The combined probability was V = 2.8 e- 04. 
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1 Introduction 

The nearby AGN Mkn 501 (z = 0.034) of the BL-Lac class is known to be one of the few TeV 7-ray 
(shortcut 7) sources. The source has been discovered in 1995 by the Whipple collaboration, Quinn et al. 1996, 
at a flux level equivalent to a few % of that of the Crab nebula and has been confirmed 1996 by the HEGRA 
collaboration, Bradbury et al. 1997. 

Early in 1997 a large increase in the TeV 7-flux has been observed by at least 7 groups, see Protheroe 1997 
for an overview. Observations covered the time from February until October. The flux was highly variable 
with a mean value of around 3 times and flares up to 10 times of that of the Crab. 

Since 1996 Mkn 501 has been monitored in the 2-12 (10) keV X-ray band by the all sky monitor (ASM) of 
RXTE (Remillard 1997). The data are available as public domain (RXTE 1999) and presented either on rates 
per few hours or smoothed. The Mkn 501 X-ray data coverage is nearly continuous (with small gaps) and 
shows also a significant increase until mid 1997 and then a falling tendency. Again, significant fluctuations 
are observed, but much less dramatic compared to the TeV flux variation. 

Already at the XXV ICRC, Durban, two groups (TA and HEGRA) reported about indications for peri- 
odicity in the TeV lightcurve, although no full analysis has been carried out. One of the problems was the 
interruption of the observations during the periods of moonshine. In the following the procedure and results 
of a QPO analysis to both the TeV and X-ray lightcurve will be presented. 

2 The data 

The TeV data were recorded with the 6 air Cherenkov telescopes of the HEGRA experiment located on 
the canary island La Palma. Four telescopes formed a stereo system (SCT) providing high statistics results 
from 110 hours of observation time. Less precise data were taken with the two standalone CTs. The latter 
instrument's data correspond to 300 hours observation, in part overlapping with the SCT data. A sizable 
fraction of the data recorded with one of the standalone telescopes, CT1, was taken during moonshine, thus 
providing flux measurements when no other telescope was operational, see Aharonian et al. 1999a, Aharonian 
et al. 1999b, Kranich et al. 1999. 

The used lightcurve above 1.5 TeV energy is shown in Fig.f^ Typical observation duration during summer 
1997 was 3-4 hours/night whenever atmospheric conditions were acceptable. Fig. ffl also shows the X-ray 
lightcurve between 2 to 10 keV from RXTE-ASM. For clarity the smoothed data are shown while in the 
subsequent analysis the unsmoothed data (~ 3000 data points) were used. 

For the analysis the data between MJD 50545 and MJD 50661 were used. As starting date we took the 
onset of the high state while the second limit was taken where the X-ray curve was interrupted for some time. 
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Figure 1: Lightcurve above 1.5 TeV energy; also shown is the X-ray lightcurve between 2 to 10 keV from 
RXTE-ASM. The one day averaged RXTE data are shown for clarity. 



Additional reason to end the analysis at MJD 50661 was the fact that the later TeV data were taken only at 
large zenith angles and rather short observation times per day. Also the source seemed to return to a lower 
state. 

3 The periodicity analysis 

As mentioned in the introduction and visible in Fig. [2 the TeV lightcurve shows indications of a periodic 
modulation. In order to test for a QPO structure we used the formalism developed by Lomb and Scargle 
(Lomb 1976, Scargle 1982) which is able to derive power spectra from unevenly sampled data. The Lomb 
normalized periodogram (nP) is a modification of the classical Periodogram (defined within the formalism of 
discrete Fourier transformations) and gives the spectral power as a function of angular frequency lo = 2tt/. 
The nP is defined as: 
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The hj are the individual flux values measured at time tj, and the sums extend over all N data points, h is the 
mean value of the hj . 

This definition of r makes Pj\r ( lo ) completely independent of a constant time shift of the data. It can be 



shown (Scargle 1982) that the normalization on the variance forces Pjy ( k> ) to have an exponential probability 
distribution if the data values are independent Gaussian random values: 



V(P N O) > z) = exp(-z) . 
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The results from the analysis, carried out independently for the TeV and X-ray data set, are shown 
in Fig. |2] Both, the TeV and X-ray data 
show a significant peak at a frequency of 
0.044 dr 1 which corresponds to a 23 day 
period. The 2nd harmonic of this pe- 
riod is visible in the X-ray data but not in 
the TeV data. Note that the high power 
value in the X-ray data at / = 1 / T = 
0.0088 d^ 1 corresponds to the overall be- 
havior of the average count-rate during 
1997 (increasing/decreasing count rates be- 
fore/after MJD ~ 50625). 
From the power value for the 23 days pe- 
riod one can calculate the corresponding 
probability against the null hypothesis of 
independent Gaussian random noise. But Figure 2: Power spectra for TeV (solid line) and x-ray data (dashed line), 
since it is known that BL Lac objects show the phenomenon of flaring it is more meaningful to test against 
the null hypothesis of independent, random distributed flares (so-called shot noise model, see deJager 1999 
and reference therein). In order to state probabilities within this framework it is necessary to renormalize nP. 
This can be done by calculating the mean Fourier power from the shot noise model. If one assumes flares of 
Gaussian shape with standard deviation a the mean Fourier power becomes: 
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Here A ex a 2 which is the mean squared amplitude of the flares. A maximum likelihood analysis of the 
derived power spectra of Fig.EJ but excluding the 23 day period and the first frequency / = 1 / T, gives the 
parameters A = 3.8, a = 0.41 d (TeV range) and A = 3.8, a = 0.35 d (X-ray range) (deJager 1999). The 
final Fourier power is then calculated as: 



P(uj) = P N (u;)/(P f (u)) 
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The distribution of P ( lo ) is shown in Fig.|3] The power values now follow the same exp ( — z ) distribution 
as in the case of independent Gaussian random noise (Equ. |3j>. Again the 23 day period points have big, but 
less significant deviations from this distribution. 

The final significance for the 23 day period is derived from the power spectrum of the combined X-ray and 
TeV data. Here we use the fact that if V = exp ( — z ), for each single frequency 

(6) 
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is chi-square distributed with 4 degrees of freedom. The resulting x 2 values together with the corresponding 
probability are shown in Fig. |3] For the 23 day period a probability V = 4.9 e-06 corresponding to 4.4 a is 
derived, which reduces to V = 2.8 e- 04 or 3.5 a after taking all 58 independent trial frequencies (between 
1/T and 0.5 per day) into account. 




Figure 3: Number of frequencies with a given 
(renormalized) power vs. power for the TeV (black 
squares) and X-ray data (open circles). Points with- 
out errors consist of one single frequency. 



Figure 4: The combined power spectrum (X-ray and 
TeV data) using proper normalization (see text). The 
probability of the 23 day period, after multiplication 
with the number of trials, corresponds to 3.5 a . 



4 Results and Conclusions 

The power spectrum shows a significant peak at 23 days in both the TeV and X-ray data samples. This 
result is in accordance with the significant correlation of 0.64 which has been observed between the TeV 
and X-ray data (Aharonian 1999b). Using the framework of the shot noise model the combined significance 
becomes 3.5 a. In order to avoid/reduce aliasing effects, introduced by the data gaps during the moon period 
in the TeV range, moon observations have proven to be of great importance. 
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